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Nomenclature 

A system dynamics matrix (n x n) 

B input distribution matrix ( n x m) 

C state-to-output distribution matrix 

(Z x n) 

C e relative modal control effectiveness 

matrix (n f x m) 

c unit-length steady-state input 

perturbation m-vector 

c Po normalized relative modal control 

factor, or element of C e (£ — 1 to 
n', j = 1 to m) 

D input-to-output distribution matrix 

(Z x m) 

i imaginary number, y/—l 

J objective function 

Z number of output perturbation 

quantities 

m number of input perturbation 

variables 

n number of state perturbation 

variables 

n f number of system modes 

O relative output significance matrix 

(Z x n f ) 

Ok i normalized relative modal observ- 

ability factor, or element of O 
(k = 1 to Z, £ = 1 to n!) 

p roll rate perturbation in airplane 

body axis, deg/sec 

Q steady-state input-to-output distri- 

bution matrix ( / x m) 

q*. m-vector formed by the fcth row of 

Q 

q pitch rate perturbation in airplane 

body axis, deg/sec 

r yaw rate perturbation in airplane 

body axis, deg/sec 

S optimization matrix (m x m) 

t time, sec 

u input or input perturbation m- 

vector 


uj input perturbation variable (j — 1 

to m) 

V nonsingular transformation matrix 

between scaled and modal coordi- 
nates (n x n) 

V t airplane total velocity perturbation, 

ft /sec 

W diagonal weighting matrix used for 

scaling 

x state perturbation n- vector 

state perturbation variable (i — 1 to 
n) 

y output perturbation /-vector 

output perturbation variable (k = 1 
to Z) 

a angle-of- attack perturbation, deg 

(3 angle-of-sideslip perturbation, deg 

r modal controllability matrix (n x m) 

7 ij element of T (i = 1 to n, j = 1 to 

m) 

6 perturbation of control surface, 

throttle lever, or thrust vector, deg 

0 Euler pitch angle perturbation, deg 

A modal dynamics matrix, block 

diagonal (n x n) 

Aj eigenvalue or element of A (i = 1 to 

n) , rad/sec 

£ state perturbation n- vector in 

modal coordinates 

real part of eigenvalue A^, rad/sec 

modal observability matrix (Z x n) 

(ffc i element of $ (k = 1 to Z, i = 1 to n) 

(j) Euler bank angle perturbation, deg 

^ Euler heading angle perturbation, 

deg 

a imaginary part of eigenvalue Aj, 
rad/sec 

Subscripts: 
a aileron 

as antisymmetric 

/ trailing-edge flap 


ill 




preceding page blank not filmed 


n 


leading-edge flap 
r rudder 

s horizontal stabilator 

ss steady state 

sy symmetric 

th engine throttle lever 

tv thrust vectored 

Superscripts: 

s scaled quantity 

T transpose of matrix [•] 

* estimated best value for scaling 

— 1 inverse of matrix [•] 

( ) derivative with respect to time, d/dt 

Abbreviations: 

D.R. Dutch roll 

Phug. phugoid 

S.P, short period 

Spir, spiral 


Abstract 

An engineering approach for analyzing airplane control and output 
characteristics is presented. State-space matrix equations describing 
the linear perturbation dynamics are transformed from physical coor- 
dinates into scaled coordinates. The equations are scaled by apply- 
ing various transformations that employ prior engineering knowledge 
of the airplane physics. Two analysis techniques are then explained. 
Modal analysis techniques calculate the influence of each system in- 
put on each fundamental mode of motion and the distribution of each 
mode among the system outputs. The optimal steady-state response 
technique computes the blending of steady-state control inputs that op- 
timize the steady-state response of selected system outputs. An exam- 
ple airplane model is analyzed to demonstrate the described engineering 
approach. 


Introduction 

For the early phases of control system design 
and analysis, engineers typically use a model of the 
airplane dynamics represented by linear state-space 
equations. The control design cycle usually begins 
with an analysis of the basic airplane without au- 
tomatic controls, followed by an engineering deter- 
mination of the appropriate control structure and 
subsequent calculation of the control system gains. 
The engineer then analyzes and evaluates the air- 
plane model including the control system. Depend- 
ing on the results of the analysis, the engineer may 
choose to adjust the design until it can be “frozen” 
and implemented in nonlinear simulation for addi- 
tional analysis. Figure 1 depicts this process. 

As airplane designers attempt to extract more 
and more performance from each new configuration, 
controls engineers are increasingly confronted with 
airplane models that have a large number of inputs 
and/or outputs (possibly redundant) and higher or- 
der dynamics. Under these conditions, it becomes 
more difficult for the engineer to directly apply tra- 
ditional analysis techniques and produce results that 
are readily understood. 

This report presents a systematic engineering ap- 
proach for analyzing linear state-space equations. 
Many of the relations presented are commonly used 
by engineers and can be obtained either directly or 
indirectly as by-products of the design process. For 
example, the engineer frequently wishes to compare 
the sizes of elements in particular matrices. Comput- 
ing the relative sizes directly assists the engineer in 
envisioning and understanding the physics of a given 
problem. A formalized approach of scaling states, 
inputs, and outputs is presented to accomplish this. 

Although some extensions to the literature are 
proposed, the primary goal of this report is not to 


present new theory, but to emphasize the engineering 
aspects of dynamics analysis for the control system 
designer. An example is used to illustrate the types of 
observations that can be made when applying these 
analysis techniques. 



* 

To nonlinear 
simulation 


Figure 1. Linear control design and analysis process. 
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Mathematical Development 

This section describes the mathematical develop- 
ment of the engineering analysis approach. In ad- 
dition to a traditional presentation of theory, en- 
gineering assumptions or implementation issues are 
also described. A diagram showing the hierarchy of 
the analysis approach is shown in figure 2. Begin- 
ning with a state-space description of a linear system, 
the system is converted from physical coordinates to 
scaled coordinates. The scaled system is used as the 
basis for all subsequent analysis. One analysis tech- 
nique transforms the scaled system to modal coor- 
dinates and then manipulates the controllability and 
observability matrices to compute quantities that in- 
dicate the relative effect of the inputs on the modes 
and the appearance of these modes in the outputs. 
The second analysis technique uses the scaled system 
to compute the “blend” of steady-state inputs that 
“optimize” one or more steady-state outputs. 



Figure 2. Hierarchy and flow diagram for engineering analysis 
approach. 

State-Space Description of Linear 
Dynamic Systems 

The motion of many dynamic systems can be 
described by a set of equations in the form 


X = fi(x,u) 1 

y = f2(x,u) J 


( 1 ) 


quantity, but here and in the subsequent develop- 
ments they are grouped for notational compactness. 

Equations (1) can be linearized about a refer- 
ence condition. This reference condition is usually 
an equilibrium point. The resulting linear equations, 
valid for small perturbations from the reference con- 
dition, are 

x = Ax + Bu 

y = Cx + Du 

where A is the system dynamic matrix, B is the in- 
put distribution matrix, C is the state- to-output dis- 
tribution matrix, and D is the input-to-output dis- 
tribution matrix. Equations (2) are called the linear 
state-space system equations, with the variables x, u, 
and y being vectors of perturbation quantities. For 
this report, the reference condition is an equilibrium 
point; therefore, A, B, C, and D are time invariant. 



Scaled Coordinates 

Interpretation of dynamic system behavior can be 
cumbersome because of the variety of units and mag- 
nitudes of the parameters used to describe dynamic 
motions. In a wide range of airplane problems scal- 
ing of variables can make the analysis more coherent. 
For example, aerodynamicists have long known the 
benefits of using nondimensional coefficients to ana- 
lyze and compare various configurations. 

Most efforts directed at solving the scaling prob- 
lem consist of enforcing a common angular measure 
(either degrees or radians) and sometimes scaling all 
velocity elements by the total vehicle velocity at trim 
(ref. 1). This subsection presents a different tech- 
nique for scaling state perturbation variables, input 
perturbation variables, and output perturbation vari- 
ables. The purpose for this scaling concept is to pro- 
vide a common basis for comparison in the analy- 
sis. As an example, this concept results in scaling 
of the inputs (controls) on the basis of each control 
effector’s maximum authority (limit). This scaling 
concept also enables engineering knowledge of the 
airplane dynamics to be captured and used in the 
analysis. 


where x is a vector of time-varying state variables, 
u is a vector of input variables, and y is a vector of 
output variables. The input variables u can be any 
combination of force or moment generators, pilot or 
command inputs, process noise or gust inputs, and 
measurement noise inputs. In the literature a sepa- 
rate term is commonly defined for each type of input 


Any system can be scaled by applying a similarity 
transformation to its state variable representation 
(ref. 2). Such a transformation does not affect the 
eigenvalues of a system. Therefore, to scale the state 
perturbation variables, define 

x = W x x 5 (3) 


2 



where x s is the scaled state perturbation vector and 
W x is an n x n diagonal matrix, 


W x = diag (xj, X 2 , x*, . x*) (4) 

where x* is the estimated best value for scaling state 
perturbation variable so that, in general 


< 1.0 


(5) 


is true for all time. A good value for x* in most cases 
would be the predicted maximum value of the state 
perturbation variable X{. This predicted maximum 
value of the state perturbation variable should be 
chosen to define the region of linearity for the model. 
The choice of x* requires engineering judgment and 
is a way that physical understanding of the system 
being studied is incorporated in the analysis. 

The input perturbation variables are scaled on the 
basis of their respective limits or authorities. Hence, 
define an m x m diagonal scaling matrix for input 
variables as 


W u = diag(w|, u* 2 , u*, u* m ) (6) 

where u*j is the maximum anticipated size for the 
input perturbation variable uj. For controls, u * 
are the control limits or authorities and are another 
example of how the physical characteristics of the 
system are included in the analysis. The largest 
possible magnitude of an input is then unity. A 
control input of ±1 is then said to be against its stop 
(limit). Therefore, the scaling of the input variables 
is applied by letting 


u = W u u 5 (7) 

Note that for many practical airplane problems, the 
control position limit is not symmetric about the trim 
control position. This point is further described later 
in the “Example Application” section. 

Finally, the output perturbation variables are 
scaled by defining 


y = w v y s (8) 

where y s is a vector containing scaled output pertur- 
bation variables and W y is a l x l diagonal matrix, 

Wj, = diag (yj, 1 / 2 ,..., y* k , V*) (9) 

where is the maximum anticipated size for the 
output perturbation variable Like x*, y £ should 


be chosen to define the region of linearity for the 
model. 

Applying the scaling operations (eqs. (3), (7), 
and (8)) yields the system equations in scaled 
coordinates as 


x s = A s x s + B 5 u 5 1 
y s = C 5 x 5 + D s u s J 


( 10 ) 


where the scaled system dynamics, input distri- 
bution, state-to-output distribution, and input-to- 
output distribution matrices are 


A s = W~ X A W x ' 
B s = W^B W u 
C^W-'CW, 
D* = W- 1 DW„ ( 


( 11 ) 


In subsequent developments, the superscript s de- 
noting the scaled quantities is omitted for notational 
compactness. 


Modal Analysis 

It is often useful to evaluate airplane linearized 
dynamics in terms of the fundamental modes of mo- 
tion. A system description in terms of its fundamen- 
tal modes is often referred to as being in modal coor- 
dinates. In this section, the scaled system (eqs. (10)) 
is first transformed into modal coordinates. Then 
the controllability and observability matrices are ma- 
nipulated to determine the relationships between the 
system modes, system inputs, and system outputs. 

The following discussion assumes that the system 
dynamics matrix A has no repeated eigenvalues — 
a reasonable assumption for airplane dynamics. An 
eigenvalue defines a system’s stability as well as nat- 
ural frequency, damping ratio, and/or time constant 
of a given system mode. For each eigenvalue, there 
exists a corresponding eigenvector. Because eigen- 
vectors define how the physical states participate in 
a given mode of motion, modes are often identified 
by examining the system eigenvectors (ref. 3). 


Modal coordinates . Define the matrix V as a 
similarity transformation assembled from the eigen- 
vectors of A as follows. For each real eigenvalue of 
A, one column of V is set equal to the correspond- 
ing eigenvector. For each complex conjugate pair of 
eigenvalues of A, one column of V is set equal to the 
real part of the corresponding eigenvectors, and the 
adjacent column is set equal to the imaginary part 
of either eigenvector (ref. 2). This definition of V re- 
sults in a matrix composed entirely of real numbers. 


3 


A transformation between the scaled coordinates 
and the modal coordinates can be defined as 


£ = v -1 x J 


( 12 ) 


corresponding to a second-order mode indicate the 
uncontrollability of that mode. 

The modal observability matrix gives the de- 
gree to which each mode is observed in each output 
variable: 


Applying this transformation to the scaled system 
(eqs. (10)) yields an expression for the system in 
scaled modal coordinates: 


where 


£ = A£ + ITi 1 
y = <££ + Du J 

A = V -1 AV 1 
r = v _1 B > 
$ = cv J 


(13) 


(14) 


Since the matrix D contains no information about 
the system states x, the transformation to modal 
coordinates does not operate on D, the input-to- 
output distribution matrix. 

The matrices A, T, and # can be used to gain 
physical insight into the system under considera- 
tion. The modal dynamics matrix A is in a block- 
diagonal form known as modified Jordan canonical 
form (ref. 4), each block containing the eigenvalue(s) 
of a system mode. First-order modes produce a sin- 
gle element on the diagonal equal to the real eigen- 
value for that mode A*. As a result of the definition 
of V, second-order modes are represented by 2 x 2 
blocks with the real, <7*, and imaginary, parts of 
the complex conjugate eigenvalues on cross diagonals 
as in the following example: 



'¥>11 

¥>12 • 

• <Pln~ 


¥> 21 

¥?22 • 

■ <P2n 

$ - 

¥>31 

¥>32 • 

• <P3n 


■ <fin 

¥>12 • ■ 

Win - 


A column of zero elements in $ corresponding to a 
first-order mode or two columns of zero elements in 
^ corresponding to a second-order mode indicate the 
unobservability of that mode. 

Engineers may examine A, T, and $ for relative 
types of comparisons as well as absolute measures 
of influence; that is, they usually compare rows or 
columns of numbers with each other — a burdensome 
process as the size of systems under consideration 
becomes large. This report proposes to precompute 
the numbers of interest and then organize them in an 
easily interpreted fashion. 

Relative control effectiveness . Assuming that 
the terms in equations (13) have been scaled, some 
physical insight may be gained from comparing the 
terms of matrix T. Engineers who use the modal con- 
trollability matrix T during airplane dynamics anal- 
ysis typically gain the most insight from comparing 
the elements in a given row with each other. This is 
nontrivial for systems with complex modal pairs. 

To simplify this process, Lallman’s relative modal 
control effectiveness matrix (ref. 5) can be computed: 


rAi 

0 


0 0 

(72 U>2 

— (72 


0 - 

0 

0 


(15) 


Lo 0 0 ... Aj 


The modal controllability matrix T gives the de- 
gree to which each control variable can excite each 
mode: 



'711 

712 ■ 

Tim 


721 

722 ■ 

72m 

r = 

731 

732 ■ 

• 73m 


-7nl 

7n2 •• 

7nrn - 


A row of zero elements in T corresponding to a 
first-order mode or two rows of zero elements in T 


Ce — [ c e^] (18) 

For each nonzero row i in T that corresponds to a 
first-order mode £ in A, 


Ce *j ~ 


h 




Et \ 


hj 


(19a) 


where 7 ij is the element in the ith row and jih 
column in the matrix I\ If i corresponds to the first 
row in T for a controllable second-order mode £ in A, 
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The ranges of the subscripts are 


column for an observable second-order mode l in A, 


i = 1, 

2, ■■ 

n 

3 = 1, 

2, .. 

. , m 

< = 1, 

2, •• 

i 

., n 


(20) 



where n is the number of system states and eigen- 
values, m is the number of system inputs, and vl is 
the number of system modes. If the system under 
study has one or more second-order modes, n f is less 
than n. If the system has an uncontrollable mode, 
the corresponding row in C e is defined to be the null 
vector. The elements in C e are the normalized rel- 
ative modal control factors. The dimension of C e is 
n' x m, and all elements in C e are nonnegative and 
less than or equal to one. 

The relative control effectiveness matrix C e is 
used by individually looking at the matrix elements 
in each row (corresponding to a mode £). A large 
value of c eej indicates that the jth control perturba- 
tion has significant influence on the ^th mode relative 
to the other control perturbations. A small value in- 
dicates a control perturbation that has little influence 
over that particular mode. Matrix elements in differ- 
ent rows of C e cannot be compared since each row 
has been normalized individually. Applications of the 
relative control effectiveness matrix C e include aid- 
ing the engineer in determining the best set of control 
effectors for a given system. 


Relative output significance . The relative 
output significance matrix O is analogous to the rel- 
ative control effectiveness matrix defined in the pre- 
vious section. This matrix can be used to determine 
how much a particular mode can be observed in the 
various outputs of the physical system under study. 
Therefore, the relative output significance matrix O 
is defined as a matrix formed by the columns of $ 
normalized for each mode £ in A; 

o = [o M ] (21) 

For each nonzero column i in $ that corresponds to 
a first-order mode l in A, 


°k,t 


m 9 i\ 



(22a) 


where is the element in the fcth row and ith 
column of the matrix $. If i corresponds to the first 


The ranges of the subscripts are 


i = 1 , 2, .. 

n ] 



k = 1, 2, .. 


► 

(23) 

e = i, 2, .. 

n J 



where l is the number of system outputs. 

If the 


system under study has an unobservable mode, the 
corresponding column in O is defined to be the null 
vector. The elements in O are the normalized relative 
modal observability factors. The dimension of O is 
l x n', and all elements in O are nonnegative and less 
than or equal to one. 

The relative output significance matrix O is used 
by individually looking at the matrix elements in 
each column (corresponding to a mode £). A large 
value of Ok £ indicates that the fcth output strongly 
observes the fth mode relative to the other outputs of 
the system. Small values of i indicate outputs that 
weakly observe a particular mode. Matrix elements 
in different columns of O cannot be compared since 
each column has been normalized individually. Ap- 
plications of the relative output significance matrix 
O include aiding the controls engineer in determin- 
ing which sensors best observe the system modes to 
be controlled. 

Optimal Steady-State Response 

The development in this section was motivated by 
a desire to provide an analysis and design technique 
analogous to the Lallman pseudo-control method 
(ref. 5). For airplane configurations with a large num- 
ber of redundant control effectors, Lallman’ s method 
can be used to allow computation of feedback gains 
for a small number of pseudo controls (ref. 6). Lall- 
man’s method computes a controls blending vector to 
maximize and/or minimize the response of a set of 
system modes. It is conceptualized for control strate- 
gies that decouple certain modes. Several airplane 
dynamics and control problems, including longitu- 
dinal rigid-body motion, require a different control 
strategy. For this class of problems, a controls blend- 
ing method is sought that would operate on the state 
perturbations of a system, rather than on the system 
modes. 

Although there are many potential solutions to 
this challenge, the author chose to develop a method 
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based on the relationship between the system input 
perturbations and the steady-state response of the 
system as measured in the output perturbations. For 
this method to have validity, some proportional rela- 
tionship must exist between the steady-state and the 
transient response of the physical states to the system 
inputs. Therefore, this method yields meaningful an- 
swers only for certain types of airplane dynamics and 
control problems. 

The technique described in this section should be 
used with caution for airplanes with unstable modes. 
This is perhaps intuitive since any unstable mode 
(and at least one of its associated physical states), 
when excited by a system input, results in an infinite 
steady-state output perturbation response. It is also 
worth noting that the “steady-state” output response 
predicted by a linear model may lie outside the 
model’s region of linearity (validity). The engineer 
must be aware of the above facts and use judgment. 

The steady-state response of the system outputs 
due to a given set of system inputs can be determined 
as follows. The existence of a steady-state condition 
implies that 

^ss = 0 (24) 

For aircraft, the steady-state condition typically 
restricts the model’s reference condition to non- 
accelerating flight. Combining the steady-state con- 
dition (eq. (24)) with system equations (10) yields 
the steady-state scaled system equations 

0 = Ax ss -j- Bu ss (25a) 

y ss = Cxgy 4“ Du ss (25b) 

where A, B, C, and D are the scaled system dynam- 
ics, input distribution, state-to-output distribution, 
and input-to-output distribution matrices defined in 
equations (11). 

The following development assumes that the sys- 
tem under study does not have any zero eigenvalues 
(singularities). When singularities exist, the system 
must be altered by deleting the physical state that 
causes A to be singular (heading angle \p is a typical 
example) . 

Having stated this caveat, equation (25a) can be 
rewritten as 

x ss = — A -1 Bu ss (26) 

Substituting equation (26) into equation (25b) yields 

Yss — CA *B 4- (27) 

The values in the scaled steady-state output pertur- 
bation vector y ss are directly related to the values 


in the scaled steady-state input perturbation vector 
Ugg. Therefore, the matrix computed in equation (27) 
is an indication of the steady-state control power of 
each system input as measured in the system outputs. 
The goal is to determine the steady-state control dis- 
tribution that maximizes a combination of selected 
outputs while minimizing a combination of undesired 
outputs. 

The distribution matrix between the steady-state 
inputs and the steady-state outputs is defined as 

Q = — CA -1 B + D (28) 

Then let 

= c u (29) 

where c is a unit-length vector of relative input 
magnitudes and u is a scalar pseudo-control variable. 
Therefore, equation (27) becomes 

yss = Q cii (30) 

Let q represent a column vector formed by the fcth 
row of Q. Form an objective function 

j = c T ( Y. q E q*tf 


Outputs 

Outputs 


to max. 

to min. 






- 5 E 

qjtq I , 

-2 
C U 

(31a) 

Outputs 

remaining 

) 




which is to be maximized irrespective of the value of 
u and subject to the constraint 

c T c = 1 (31b) 

The constraint in equation (31b) is chosen so that 
for the scaled system (eqs. (10)), the controls have 
maximum deflections of ±1. The first term of J 
is a summation of the scaled steady-state output 
perturbations to be maximized by selection of c, and 
the second term of J is a summation of the scaled 
steady-state outputs to be minimized by selection of 
c. Outputs included in the remainder term of J have 
some nonoptimal value. 

The implementation indicated by equations (31) 
was developed to allow selection of the minimum 
set of outputs to be optimized without causing nu- 
merical problems when using eigensystem decompo- 
sition software. However, this implementation is not 
without potential problems. The value chosen for 
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the remainder weighting (10 -5 ) is based on the au- 
thor’s experience with aircraft models using the scal- 
ing techniques described in this paper. This weight- 
ing must be sufficiently large to prevent numerical 
problems associated with the computer algorithm be- 
ing used. Conversely, if the weighting is not suf- 
ficiently small relative to the elements in the row 
vectors q*., the “optimization” may correspond to 
one of the remainder outputs — an undesirable result. 
Proper scaling should prevent this problem, but the 
reader is advised nevertheless. 

An equivalent formulation to equations (31) can 
be written as 


If this is not satisfactory, the above development can 
be extended in a straightforward manner by applying 
different weights to distinguish between the steady- 
state outputs. 

Display of Computed Functions 

The previous section describes the computation 
of various vectors and matrices which the control 
system engineer can use to gain physical insight 
into the dynamic system under study. Graphic dis- 
play of these quantities can further aid the engineer 
in quickly analyzing and understanding a dynamic 
system. 


m “(^) =max (w) 


(32) 


where 


S 


£ ^ 1 - £ q * q f 


Therefore, a format was developed for graphically 
presenting the various vector and matrix quantities 
which are computed using the analysis techniques 
described previously. Both vectors and matrices are 
presented in a bar chart format as shown in figure 3. 
For matrices, multiple bar charts are plotted for the 
set of vectors that define a given matrix. 


Outputs 
< to max. 


Outputs 
to min. 


\ 


+ 10 5 £ q fc q[ 


Outputs 

remaining 


/ 


(33) 


for which the solution is well-known (ref. 7) and 
described by 


max 


'c^Sc 

c^c 


— A max (S) 


(34) 



Therefore, the solution for c is the eigenvector which 
corresponds to A ma x(S), the maximum eigenvalue of 
the matrix S. 

In the development described in equations (30) 
through (34), it has been implicitly assumed that 
the steady-state system output perturbations to be 
maximized or minimized were of equal importance. 


Figure 3. Bar chart format. 

The bar chart vertical axis is the value of the 
vector elements. The horizontal axis indicates the 
elements of the vector. Each bar is located over an 
integer tick mark indicating a vector element. The 
height of the bar as measured against the vertical 
axis denotes the value of the vector element. 
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Example Application 

In this section, an example application is described to illustrate some of the observations that can be made 
using the analysis techniques given in this report. 

The example is a modern high performance airplane with a large number of redundant control effectors 
and hypothetical thrust vectoring in trimmed level flight at an angle of attack of 20°. Figure 4 is a drawing of 
the example high performance airplane. A thrust-vectoring system consisting of turning vanes configured as 
shown in figure 5 is postulated. For each engine, two turning vanes are used to deflect the exhaust flow. Each 
set of turning vanes is oriented so that each upper vane is rotated upward 48° from the airplane’s horizontal 
axis. This thrust-vectoring system can generate control moments in the pitch, roll, and yaw axes. Note that 
yaw control moments cannot be generated without also generating roll control moments, and vice versa. 

The analysis techniques were implemented using the MATRIX x software package. (MATRIX x is a registered 
trademark of Integrated Systems, Incorporated, Santa Clara, California.) The MATRIX x programming 
capability was used as the language for implementing the modal analysis techniques. Reference 8 lists some of 
the software used to analyze the examples. 

The following linear system perturbation equations are for the example airplane at a trimmed, straight and 
level flight condition: altitude of 5000 ft, angle of attack of 20° and total velocity of 235.1 ft /sec. The output 
perturbation variables are postulated by the author. 

State-Space Description 

The linear system perturbation equations for the example airplane are given in equations (2) where 


A = 


B = 


C = 


X 

= {Vt 

a q 

6 (3 p 

r 4)} 

T 






(35a) 

y 

= {V t 

a 0 

p q r 

<P 0} 

T 






(35b) 

u 

= { 6s *y 

^Tl S y 

6 fsy Ky 

<5tV S y 

% 

^ s as fyas ^as ^Vas 

1 T 
6 r } 


(35c) 

- 0.069 

- 0.360 

0 

- 0.561 

0 

0 

0 

0 ' 





- 0.061 

- 0.303 

0.986 

0 

0 

0 

0 

0 





- 0.006 

- 0.407 

- 0.236 

0 

0 

0 

0 

0 





0 

0 

1.0 

0 

0 

0 

0 

0 




(35d) 

0 

0 

0 

0 

- 0.091 

0.342 

- 0.940 

0.129 




0 

0 

0 

0 

- 7.89 

- 0.705 

0.669 

0 





0 

0 

0 

0 

- 0.130 

0.013 

- 0.088 

0 





0 

0 

0 

0 

0 

1.0 

0.364 

0 





- 0.120 

0.104 

- 0.095 

0 

- 0.026 

0.168 

0 

0 

0 

0 

0 ' 


- 0.070 

- 0.006 

- 0.047 

- 0.025 - 

0.016 

- 0.015 

0 

0 

0 

0 

0 


- 1.69 

- 0.091 

0.145 

- 0.133 - 

0.529 

0.028 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(35e) 

0 

0 

0 

0 

0 

0 

- 0.003 

0 

- 0.002 

- 0.015 

0.014 

0 

0 

0 

0 

0 

0 

2.95 

3.48 

2.52 

0.223 

0.382 


0 

0 

0 

0 

0 

0 

0.005 

- 0.061 

- 0.095 

0.420 

- 0.290 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



'1 0 0 0 0 0 0 0 ' 

0 1 0 0 0 0 0 0 

0 0 0 0 1 0 0 0 

0 0 0 0 0 1 0 0 

0 0 1 0 0 0 0 0 

0 0 0 0 0 0 1 0 

0 0 0 0 0 0 0 1 

.0 0 0 1 0 0 0 0 . 


D = O 


(35f) 
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All input perturbations (u) are in degrees. The total velocity ( Vt ) is in ft/sec; the velocity vector angles 
and Euler angles (a, /?, <f > , and 9) are in degrees; the body axis rates (p, g, and r) are in deg/sec. 

Scaled Coordinates 

The diagonal weighting matrices defined in equations (4), (6), and (9) must be determined to apply the 
scaling operations of equations (3), (7), and (8). The following guidelines are used to select the maximum 
expected values of the state perturbation variables used in equation (4): 

1. Define V t * to be 10 percent of V^, the total airplane velocity at trim. 

2. Define a* and (3* to maintain the airplane in the region of approximate linear aerodynamics about trim. 

3. Define 6* to be 5° and <f>* to be 10° (linear region for Euler angles about trim.) 

4. Define p* to equal <j > * in equivalent units (to maintain linearity in Euler angles for 1 sec). 

5. Define q* and r* to equal one- fourth of p*. 

These guidelines result in 


W, = diag(v;, a*, g*, 9\ /T, p\ r\ 0*) 

= diag(24, 5, 2.5, 5, 2.5, 10, 2.5, 10) (36) 

These same guidelines are used to determine the maximum expected values of the output perturbation variables 
used in equation (9): 


W y = diag(V?, p\ r\ 0*, 0*) 

= diag(24, 5, 2.5, 10, 2.5, 2.5, 10, 5) (37) 

Table 1 lists the guidelines used to select the maximum anticipated values of the input variables used in 
equation (6). Several issues associated with selection of the maximum anticipated size for airplane control 
input variables require engineering judgment. Each engineer may apply this judgment differently. During 
the analysis, the author developed two “rules” to describe how engineering judgment influenced the scaling 
methodology. The author offers these rules with the knowledge that they may be the subject of future debate. 

Table 1. Guidelines for Scaling Input Variables 


Variable 

Description 

Lower limit, 
deg 

Upper limit, 
deg 

Trim value, 
deg 

Maximum 
anticipated 
size, deg 

£s 8y 

Sym. horiz, stabilator 

-24 

10.5 

-6.5 

6.9 


Antisym. horiz. stabilator 

-10 

10 

0 

10 

finny 

Sym. leading-edge flap 

-3 

34 

26.6 

7.4 

fife ty 

Sym. trailing-edge flap 

-8 

45 

7.8 

15.8 

fifes 

Antisym. trailing-edge flap 

-8 

8 

0 

8 

fi(Lny 

Sym. aileron 

-25 

42 

0 

1 

fi&e& 

Antisym. aileron 

-25 

25 

0 

24 

£tV 8 y 

Sym. thrust vectoring 

-25 

25 

0 

12.5 

^tVas 

Antisym. thrust vectoring 

-25 

25 

0 

12.5 

fitb. 

Throttle lever 

31 

106.5 

89.6 

16.9 

fir 

Rudder 

-30 

30 

0 

30 
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First, most of the symmetric and antisymmetric controls cannot be at their respective limits at the same 
time. As an example, antisymmetric thrust- vectoring vanes cannot be deflected when a full authority (25°), 
symmetric thrust-vectoring input is made. Therefore, the maximum anticipated values for these control 
effectors are reduced from the physical limits to allow simultaneous full authority deflections of the symmetric 
and antisymmetric controls. For some control effectors, such as the thrust-vectoring vanes, the maximum 
anticipated values allow equal authorities for the symmetric and antisymmetric controls. For other control 
effectors, such as the aileron, the maximum anticipated values favor one control (in this case, antisymmetric) 
over another. 

Second, the most conservative maximum anticipated value is used for nonsymmetric control authorities 
about trim. An example of this rule, combined with the first rule, is the horizontal stabilator. The -6.5° trim 
position for symmetric horizontal stabilator 6 Say permits a deflection from trim of —17.5 and +16.9 . Applying 
the first rule reduces the permissible deflection to -7.5° and +6.9°. Therefore, the maximum anticipated value 
for 6 Ssy is 6.9°. 

These guidelines for scaling the input variables result in 


W u = diag (6* 6, 


Ss y * Tlsy : 


n sy- 


ttgy > 


y tv. 


sy ■ 


^th’ ^as’ 


^aas’ 


u tVas> 


6 * r ) 


diag(6.9, 7.4, 15.8, 1, 12.5, 16.9, 10, 8, 24, 12.5, 30) 


(38) 


Modal Analysis 

Modal coordinates. The system (eqs. (35)) is transformed into modal coordinates with the following 
result: 


r = 



■ 0 

- 0.6543 

0.0776 

0 

0.1599 

- 0.1046 

0 

0 _ 



0 

0.1573 

0.0251 

0 

0.0803 

- 0.6561 

0 

0 



0 

- 0.3048 

0.1155 

0 

0.8983 

0.0305 

0 

0 



0 

0.3692 

0.8886 

0 

- 0.2764 - 0.6082 

0 

0 


v = 

- 0.0705 

0 

0 

0.0698 

0 

0 

0.6268 

0.3682 



0.1675 

0 

0 

- 0.2886 

0 

0 

- 0.6451 

0.5948 



- 0.2774 

0 

0 

0.1141 

0 

0 

- 0.0055 

0.0737 



. - 0.9434 

0 

0 

0.9481 

0 

0 

0.4370 

0 . 3536 . 



' - 0.1507 

0 

0 

0 

0 

0 

0 

0 

1 


0 

- 0.0053 

0.1693 

0 

0 

0 

0 

0 


0 

- 0.1693 

- 0,0053 

0 

0 

0 

0 

0 


0 

0 

0 

- 0.2934 

0 

0 

0 

0 

A = 

0 

0 

0 

0 

- 0.2989 

0.6026 

0 

0 


0 

0 

0 

0 

- 0.6026 - 0.2989 

0 

0 


0 

0 

0 

0 

0 

0 

- 0.2196 

1.554 


0 

0 

0 

0 

0 

0 

- 1.554 

- 0 . 2196 . 


0 

0 

0 

0 

0 

0 

1.0968 

2.400 

8.333 

- 13.30 

22.72 

- 1.324 

- 0.1306 

0.3395 

- 0.0164 

- 0.7594 

- 0.1381 

0 

0 

0 

0 

0 

- 1.735 

- 0.0917 

0.4984 

- 0.0158 

- 1.001 

0.1480 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.7685 

2.054 

7.538 

- 13.02 

22.03 

- 5.378 

- 0.3311 

1.053 

- 0.0625 

- 3.056 

0.1449 

0 

0 

0 

0 

0 

- 0.8951 

- 0.0628 

0.4556 

- 0.0047 

- 0.5342 

0.0665 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

- 1.790 

- 1.769 

- 4.083 

0.6504 

- 2.002 

0 

0 

0 

0 

0 

0 

3.078 

3.081 

7.070 

- 1.395 

4.046 


(39a) 


(39b) 


(39c) 
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$ = 


(39d) 


■ 0 

-0.6543 

0.0776 

0 

0.1599 

-0.1046 

0 

0 

0 

0.1573 

0.0251 

0 

0.0803 

-0.6561 

0 

0 

-0.0705 

0 

0 

0.0698 

0 

0 

0.6268 

0.3682 

0.1675 

0 

0 

-0.2886 

0 

0 

-0.6451 

0.5948 

0 

-0.3048 

0.1155 

0 

0.8983 

0.0305 

0 

0 

-0.2774 

0 

0 

0.1141 

0 

0 

-0.0055 

0.0737 

-0.9434 

0 

0 

0.9481 

0 

0 

0.4370 

0.3536 

_ 0 

0.3692 

0.8886 

0 

-0.2764 

-0.6082 

0 

0 


The eigenvalues and eigenvectors of the system in equations (39) are listed in tables 2 and 3. All the system 
modes are stable. Note that the phugoid and short period modes are closer in natural frequency than is often 
the case for airplanes. The roll and spiral modes have similar characteristics because of the relatively high trim 
angle of attack. 


Table 2. Eigenvalues of System 


Eigenvalues 

Natural 

frequency, 

sec" 1 

Damping 

ratio, 

rad/sec 

Time 

constant, 

sec 

Period, 

sec 

Time-to 

half-amplitude, 

sec 

Description 

o 

tO 

d 

1 

II 



6.64 


4.60 

Spiral 

A 2 = -0.0053 ±i0. 1693 

0.169 

0.031 


37.1 

131 

Phugoid 

x 3 = -0.2934 



3.41 


2.36 

Roll 

A 4 = -0.2989 ± i0.6026 

0.673 

0.444 


9.34 

2.32 

Short period 

A 5 = -0.2196 ± i 1 .554 

1 1.57 

0.140 


4.00 

3.16 

Dutch roll 


Table 3. Eigenvectors of System 


State 

Eigenvectors | 

Spiral 

Phugoid 

Roll 

Short period 

Dutch roll 

Vt 


0 'j 


(-0.65±i0.08) 


0 ] 


0.16 ±<0.10 ] 


0 ] 


a 


0 


0.16 ± <0.03 


0 


0.08 ± <0.66 


0 


Q 


0 


-0.30 ±<0.12 


0 


0.90 ± <0.03 


0 


9 


0 


0.37 ± t0.89 


0 


-0.28 ± <0.61 


0 


P 


-0.07 

y < 

0 

y < 

0.07 

> < 

0 

► i 

0.63 ± i0.37 


V 


0.17 


0 


-0.29 


0 


-0.65 ± <0.59 


T 


-0.28 


0 


0.11 


0 


-0.01 ±<0.07 




k -0.94 , 


0 


, 0.95 i 


, 0 


k 0.44 ± <0.35 j 



Relative control effectiveness . The relative control effectiveness matrix C e defined in equations (19) 
is a collection of normalized relative modal control factors which enable comparison of the system inputs’ 
effectiveness for a given mode: 


^sy 

bn sy 

bfsy 

ba sy 

^tVgy 

<5th 

^*as 

6 f 

Jas 

ba as 

^tVas 

6r 



0 

0 

0 

0 

0 

0 

0.040 

0.087 

0.300 

0.480 

0.819 

Spir. 


0.839 

0.061 

0.232 

0.009 

0.483 

0.078 

0 

0 

0 

0 

0 

Phug. - 


0 

0 

0 

0 

0 

0 

0.029 

0.077 

0.282 

0.486 

0.823 

Roll , 

bio not 

0.853 

0.053 

0.180 

0.010 

0.486 

0.025 

0 

0 

0 

0 

0 

S.P. J 

, compare 

0 

0 

0 

0 

0 

0 

0.333 

0.332 

0.762 

0.144 

0.422 

D.R. 



Compare (40) 
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Figure 6 shows plots of the row vectors of C e for the example airplane, the values for which are indicated m 
equation (40). The first six columns of matrix C e correspond to the longitudinal (symmetric) inputs, whi e 
the remaining columns correspond to the lateral-directional (antisymmetric) inputs. 

Recall that the rows of this matrix correspond to the system modes and the columns correspond to the 
system inputs. For example, the lower right element in C e (0.422) indicates the effectiveness of the rudder 
on the Dutch roll mode relative to the other controls. Also recall that comparisons cannot be made between 
matrix elements in different rows of C e . For example, one cannot draw the conclusion that the rudder b r is 
more effective over the roll mode (0.823) than the Dutch roll mode (0.422). 

Of all the controls, the symmetric stabilator 8 Ssy is the most effective in exciting both the phugoid and the 
short period mode; this is observed by comparing values in the second and fourth rows, respectively, of C e . 
The rudder 6 r is the most effective control in exciting the roll mode; this is observed by comparing values in 
the third row of C e . The rudder is less effective in exciting the Dutch roll mode than antisymmetric aileron 
6 ai<s ; this is seen by comparing values in the fifth row of C e - 


Relative output significance. The relative output significance matrix O defined in equations (22) is a 
collection of normalized relative modal observability factors which enable comparison of the system outputs 
ability to observe the motion of a given mode: 


Spir. 

Phug. 

Roll 

S.P. 

D.R. 



0 

0.540 

0 

0.145 

0 

v t 


0 

0.130 

0 

0.503 

0 

a 


0.071 

0 

0.070 

0 

0.571 

13 i 


0.168 

0 

0.289 

0 

0.689 

p 

Compare 

0 

0.267 

0 

0.684 

0 

q 

0.277 

0 

0.114 

0 

0.058 

r 


0.943 

0 

0.948 

0 

0.442 

<S> 


0 

0.788 0 0.508 

i ► 

Do not compare 

0 

9 



(41) 


Figure 7 is a plot of O for the example airplane, the values for which are indicated in equation (41). The 
first, second, fifth, and eighth rows (vertical bars in fig. 7) correspond to the longitudinal outputs, while the 
remaining rows (vertical bars in fig. 7) correspond to the lateral-directional outputs. 

Recall that the rows of this matrix correspond to the system outputs and the columns correspond to the 
system modes. For example, the third element in the rightmost column of O (0.571) indicates the appearance 
of the Dutch roll mode in the angle-of-sideslip output (3 relative to this mode’s appearance in the other system 
outputs. Also recall that comparisons cannot be made between matrix elements in different columns of O. For 
example, one cannot draw the conclusion that angle of attack a better observes the short period mode (0.503) 
than the phugoid mode (0.130). 

The phugoid mode is primarily observed in 9 and Vf, this can be seen by comparing values in the second 
column of O. The short period mode is mostly observed in q, 9 , and a; this is seen by comparing the elements 
of the fourth column of O. Of all the outputs, the Euler roll angle <j> best observes the roll mode; this is seen 
by comparing values in the third column of O. The roll angle is less effective at observing the Dutch roll mode 
than both the roll rate p and the sideslip angle f3 ; this is seen by comparing values in the fifth column of O. 
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Figure 6. Relative control effectiveness matrix C e for example airplane. 





Value Value Value 


Spiral mode 



1 . 0 - 
. 8 - 
» . 6 - 
> . 4 - 


Phugold mode 


. 2 - 

0- 


Vt a p P 


<$> e 


Roll mode 



Short period mode 

1.0 1 

. 8 - I 

I- 6 - n 


v t a P P q r <t> 6 


Dutch roll mode 



Figure 7. Relative output significance matrix O for example airplane. 





Optimal Steady-State Response 

The system’s steady-state input-to-output distribution matrix Q is 


°s B y 

°”sy 

d fs y 

Odsy 

Otv sy 

6.00 

0.338 

-1.79 

0.054 

3.46 

-6.15 

—0.356 

1.25 

-0.069 

-3.50 

0 

0 

.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.100 

0.301 

-0.279 

0.012 

0.078 


^th 

^Sas 

*/. 


^tvag 

S r 

-0.444 

0 

0 

0 

0 

0 

0.267 

0 

0 

0 

0 

0 

0 

1.30 

0.997 

1.64 

2.41 

-3.30 

0 

0.148 

0.321 

1.11 

-1.76 

3.02 

0 

0 

0 

0 

0 

0 

0 

-1.63 

-3.52 

-12.2 

19.4 

-33.1 

0 

-3.11 

-7.11 

-24.9 

40.6 

-69.4 

1.10 

0 

0 

0 

0 

0 


Vt 

a 

0 

P 

Q 

r 

4 > 

6 


( 42 ) 


where the matrix Q is defined in equation (28). The rows in Q correspond to the system outputs and the 
columns correspond to the system control inputs. Note that the row in Q corresponding to the steady-state 
pitch rate perturbation output is a null vector. This is a direct result of the model structure and the 
imposition of the steady-state condition. Since the model was trimmed in straight and level flight 6 = a. Since 
^ss — 0, then < 7 SS = 0 for this model. 


The blending of input perturbations to achieve an optimal steady-state response for a thrust-vectored fighter 
airplane can be computed as indicated in equations (30) through (34). Figure 8 plots the scaled (white bars) 
input control deflections that maximize response in the output perturbation, but do not minimize any 

steady-state outputs. The result is that the dominant controls are symmetric thrust vectoring and symmetric 
stabilator. 


As discussed previously, when one desires to compute the relative sizes of input perturbations using these 
analysis techniques, proper scaling of the input variables is important. The shaded bars in figure 8 show the 
result of not applying the input scaling of equation (38) during the analysis. To arrive at this result, the 
blending vector c was computed for a version of the example system model with unsealed inputs. This c’ was 
then transformed to scaled coordinates and normalized to unit length for comparison purposes. Recall that the 
white bars indicate the values computed when the system inputs are scaled. Not scaling the inputs would cause 
inappropriate utilization of symmetric thrust vectoring, trailing-edge flaps, and aileron, resulting in suboptimal 
ass response. 


1.0 



Symmetric inputs 


Figure 8. Input distribution that maximizes response (with and without input scaling). 
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Review 

The most effective longitudinal control for the 
airplane model is the symmetric horizontal stabilator 
6 Ssy as seen from the relative control effectiveness 
matrix C e . Similarly, examination of the relative 
control effectiveness matrix C e indicates that the 
most effective lateral control is the rudder 6 r and the 
most effective directional control is antisymmetric 
aileron 6 aas . 

With the exception of the spiral mode, the air- 
plane modes are observed in their conventional out- 
puts (e.g., the phugoid mode is observed in Euler 
pitch angle 6 and total velocity V£). The spiral mode 
was observed mostly in Euler bank angle <p since Eu- 
ler heading angle ip was not included in the model. 
The relative output significance matrix O indicates 
that all airplane modes are being observed; if control 
of the classical spiral mode is important, it may be 
necessary to add ip as a system output to distinguish 
this airplane model’s spiral mode from its roll mode 
for control design purposes. 

The control deflections that optimized the re- 
sponse were consistent with the control effectiveness 
computations. For this example, input scaling af- 
fected the control deflection magnitudes for optimal 
a ss response — especially for symmetric aileron 6 asy 
and symmetric thrust vectoring S{ Vsy . 

Concluding Remarks 

A systematic methodology has been outlined that 
transforms linear state-space matrix equations in 
physical coordinates to a set of matrix equations in 
scaled modal coordinates. The scaling processes are 
practical suggestions for handling the system units 
and control input limits that require an understand- 
ing of the system’s physics. Guidelines for perform- 
ing the scaling are indicated for the example airplane. 
Once the system was transformed to scaled modal co- 
ordinates, techniques were defined to analyze the ef- 
fects of the system perturbation inputs on the funda- 
mental modes of motion and the appearance of these 
modes in the system perturbation outputs. Also de- 
fined was a technique to compute a steady-state input 
perturbation vector that optimizes the steady-state 
response of selected output perturbations. Graphic 
formats were developed to display the vectors and 
matrices computed during the analysis. This was 
done to promote quick understanding of the analysis 
results for the dynamic system under study. 

These analysis techniques have been applied to an 
example thrust-vectoring airplane model. This ex- 
ample is representative of the types of dynamics that 
are encountered in current vehicle designs. The re- 


sults of this example application illustrate the under- 
standing which can be acquired using the described 
modal analysis techniques. 

The results of the analysis shown in this report 
can also be utilized directly in the control design pro- 
cess if desired. For example, the engineer can use 
the relative control effectiveness matrix to decide how 
the different control effectors should be allocated for 
various control tasks. Similarly, the relative output 
significance matrix can be used by the engineer to 
determine the appropriate set of sensors to use in ob- 
serving the system modes. For many airplane control 
system designs, both these tasks can be easily accom- 
plished using well-known heuristics instead of the ap- 
proach described in this report. However, for many 
modern airplane designs, with a large number of con- 
trol effectors (often redundant) and where active con- 
trol of structural vibration modes must be accom- 
plished simultaneously with control of the rigid-body 
modes, the relative control effectiveness and relative 
output significance analyses are of value. Finally, 
for the special case of an airplane configuration with 
a large number of redundant control effectors, the 
blending computation from the optimal steady-state 
response technique or the pseudo-control method of 
Lallman can be used to reduce the control law de- 
sign by allowing computation of feedback gains for a 
small number of pseudo controls. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
February 6, 1991 
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